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Photoassociation (PA) of ultracold metastable helium to the 2s2p manifold is theoretically inves- 
tigated using a non-perturbative close-coupled treatment in which the laser coupling is evaluated 
without assuming the dipole approximation. The results are compared with our previous study 
[Cocks and Whittingham, Phys. Rev. A 80, 023417 (2009)] that makes use of the dipole approxi- 
mation. The approximation is found to strongly affect the PA spectra because the photoassociated 
levels are weakly bound, and a similar impact is predicted to occur in other systems of a weakly 
bound nature. The inclusion or not of the approximation does not affect the resonance positions 
or widths, however significant differences are observed in the background of the spectra and the 
maximum laser intensity at which resonances are discernable. Couplings not satisfying the dipole 
selection rule | J — 1| <J'< |J + 1| do not lead to observable resonances. 

PACS numbers: 32.70.Jz, 34.50.Cx, 34.50.Rk, 34.20.Cf 
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I. INTRODUCTION 

Photoassociation (PA) is a powerful technique ex- 
ploited by researchers to probe the fundamental inter- 
actions of ultracold quantum gases 0,(2]. Of particular 
interest is PA in metastable helium where release of the 
large internal energy during collisions provides unique 
experimental investigation strategies. We have recently 
analysed the PA process [|| and obtained detailed infor- 
mation of the laser intensity dependence of the line shifts 
and widths of the resonance peaks in the PA spectra of 
ultracold metastable helium excited to the J — 1,0+ rovi- 
brational states in the 2 3 Si + 2 3 Po asymptote. Two 
variants of a full nonperturbative multichannel, close- 
coupled treatment were used to obtain the required scat- 
tering matrix S in the presence of the non-vanishing 
asymptotic radiative coupling, one based upon dressed 
states, and the other on a modified radiative coupling 
that vanishes asymptotically. Although both methods 
gave nearly identical results for the line shifts and widths, 
these peaks were superimposed on very significant back- 
grounds that differed, especially at higher laser intensi- 
ties. These significant backgrounds are a direct conse- 
quence of the weakly bound nature of the excited levels 
in metastable helium and are not found in PA studies of 
other systems. 

In common with all previous investigations of PA in 
ultracold gases, we assumed the dipole approximation 
in evaluating the matrix elements of the laser coupling. 
However, in light of the sensitivity of the background ra- 
diation loss to the form of the radiative coupling at large 
distances, we revisit the validity of this approximation 
and its applicability to weakly bound levels in general. 

The laser coupling term is given by 
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sition Ti of the ith electron is 

Here (a^) are the creation (annihilation) operators for 
a photon of angular frequency u>£ and polarization and 
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where pi 



-iKV r , and the vector potential at the po- 



where k is the wave- vector of the laser field, e and m are 
the electron charge and mass respectively and V is the 
normalization volume. The dipole approximation makes 
the assumption that exp(ife-rj) « 1. If the laser coupling 
is only significant in the region of the excited rovibra- 
tional states this is a reasonable assumption to make, as 
the outer turning points of the ultra long-range J = 1,0+ 
vibrational states of this investigation place an upper 
bound of r, < R/2 < 235 ao, where R is the interatomic 
distance. Since k = 1/3258.17 a^ 1 for the 2s 3 S - 2p 3 P 
transition, this gives 0.92 < | exp(zfe • rj)| < 1.08. 

The multichannel calculation of our previous investiga- 
tion, however, involves open channels of the metastable 
basis coupled to the excited state, and the S'-matrix 
elements and loss profiles were obtained by matching 
the asymptotic forms of the open channels at distances 
R > 10 5 ao, much greater than the interatomic ranges of 
the vibrational states. The coupling to the excited state, 
even at very large ranges where the uncoupled bound 
wavefunction was negligible, still influenced the calcula- 
tion due to the closeness of the bound levels to the excited 
state dissociation limit. This influence can be suppressed 
and the background to the PA spectra removed by artifi- 
cially deepening the potential well. The presence of this 
background introduced new features such as the elimi- 
nation of resonance peaks by saturation of the laser in- 
tensity, and it is evident that the laser coupling must be 
treated more carefully in non-perturbative calculations 
of PA profiles. 
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The validity of the dipole approximation at large inter- 
atomic separations of O(10 5 ) a is questionable, at least 
when weakly bound levels are considered. In this pa- 
per we study the effects of the non-dipole contributions 
to the PA spectra in metastable helium by firstly deriv- 
ing an expression for the matrix elements of the laser 
coupling which fully includes the exp(ife • rA factors and 
then re-do the multichannel calculations of |3J using these 
matrix elements. 

It is worth noting that the exponential term that is 
ignored in the dipole approximation has the effect of in- 
troducing a momentum transfer from the absorbed pho- 
ton to the molecule. For ultracold atoms this momentum 
transfer can be much larger than the initial momentum 
of the colliding atoms. 



II. EXACT LASER COUPLING 
A. General form 

The close-coupled equations for the two colliding atoms 
in the applied laser field involve the matrix elements of 
the laser coupling operator, H lnt between basis states of 
the general form 



ir 1 G 9 (i2)| 5 >|n > w > e A ) 



and 
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Here the molecular basis states |a), where a = {g, e}, are 
assumed to have no dependence upon the interatomic 
distance R but may still depend upon the molecular ori- 
entation in the space-fixed frame, specified by the angles 
(0,4>), and the electronic coordinates n. The laser field 
states arc denoted by \n,uj,e\), representing n photons 
of frequency u> and polarization e\, where A = 0,±1 for 
7r, polarization. We assume the laser is directed along 
the space-fixed Oz-axis. 

We consider the matrix element 



V. 



<'<j 



(e\(n - l,u,ex\H int \g)\n,(jj,e\). (6) 



As only the term involving the annihilation operator 
contributes, we have 



where 



l,u>,e\\a(\n,u),ex) = 5\^\/n 



(8) 



has been used. 

Evaluation of (0 requires the choice of a coordinate 
system. We specify the origin of the coordinate system 
to be the center of mass of the dimer (assumed to consist 



of identical nuclei), such that the position of nucleus A 
is given by —R/2 and nucleus B by R/2. The electron 
coordinates r, can then be expressed as r< = fjiR/2 + r' i 
where the operator rji has eigenvalues of =pl if electron 
i is centered upon atom A or B respectively, and r[ is 
the position vector of electron i with origin at its respec- 
tive atom. Similarly, the momentum operators become 
p{r i ) = f, i p{R)/2+p{r'A. 

The matrix element (JT)) then becomes 
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This expression can be simplified by separating the basis 
states into a rotational part which depends only upon 
(9, <j)) and the electronic part that is independent of R: 



(10) 



After action of the momentum operators, the matrix el- 
ement (e| . . . \g) inside the summation becomes 



rot I k ' £ A fjjik-R, 



+ W>nf e^ k - R l 2 {4>tW k <\^f)P(R) ■ e x W g ot ) 

+ (^\e^ ik - R ^^f\i(k ■ e x y h <\^W g ot ) 

+ (^ ot |^ 4fc - R/2 (V'e 1 |e lfc -'-''p(rO|V^ 1 ) -eA|^ ot )- ( U ) 

The first and third terms are zero since k - e x =0. In the 
second and fourth terms we can assume exp(ifc • r,-) pa 1 
as the inner product . . . \ipf) is non-negligible with 
respect to r[ only in the regions of the atomic electrons. 
As such, the second term is also zero because the ground 
and excited states are orthogonal. 

To proceed, we use the expansion for a plane wave in 
terms of Legendre polynomials 

e ^k-R/2 = J2i p (2p+l)j P (kR/2)P p (f] l cos9) 
v 

= ^P(2p + l)j p (kR/2)(fj i Y'D& (12) 



where D 3 m , m = D J m , m ((j),9,0) is a Wigner rotation ma- 
trix 0. Recalling the equality ep{r'A/m = i[H mo \ 7 d l ]/Ti, 
where d % = er^ is the atomic dipole operator for electron 
i, the matrix element becomes 



V 11 



2huje V 



^HP +1 (2p+l)j p (kR/2)D p * 



nhuj 
2^V 



J2(€ Ot \i P+1 (2p+l)3 P (kR/2)D p * 



(13) 
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In obtaining (|T3|) wc have assumed the difference between 
electronic energies E e — E g is approximately Tilo. 

The matrix element of d 1 is most easily evaluated in the 
molecular frame using spherical tensors. The expansion 
required is 



d'-e x = £(-l)"(e A )_„3 





(14) 



where the subscripts and j3 denote the spherical ten- 
sor components in the space- and molecular-fixed frames 
respectively and (e^)-^ = S\.fi- This expansion allows 
the complete separation of the rotational and electronic 
parts of the matrix element: 



ymt = 

eg 
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2^V 
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Here the two rotation matrices, D\*pD P Q, have been com- 
bined using standard angular momentum theory Q. 



Here the A and B indicate the atom to which the orbital 
configuration 7, belongs, N lt ~ l2 is a normalization factor 
and pls = wi +W2 + Li + L2 — L + Si + S% — S + N + w 
Q, where Wi refers to the atomic inversion symmetry of 
7i. The subscript '— ' indicates that the state has been 
antisymmctrized with respect to electron permutation: 



\^Lsn L n s )- 



1 

71 



[\^Lsn L n s ; ri ,r 2 ) 



\ 1 ^LSn L n s ;r 2 ,r 1 )] .(18) 



When Q L + n s = 0, the states \jfjfLSn L Q s )- must 
also be properly symmetrized with respect to a v , the re- 
flection operator of the electronic wave function through 
a plane containing the intermolecular axis, however this 
docs not affect the action of (f)i) p< ^a an d wm be ignored. 
For the 2s2s and 2s2p states, the symmetrization reduces 
to 

\2s2s) = X - [1 + (-I)— 5 ] |(2 S ) A (2 S ) B 0S0ft s )_ (19) 



and 



\2s2 V ) = -^=[\(2 S ) A (2 P ) B isn L n s ) 



+ (-i) 1 - s+w \(2 P r(2srisn L n s )-\ (20) 



B. Explicit form 

The appropriate explicit basis states are the hybrid 
Hund case (c) states 



|a) 
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where L a , S a and j a = L a + S a are the orbital, spin 
and total angular momenta of atom a and 7 = {7172} 
where j a = {^( a L a S a } with all other necessary atomic 
quantum numbers specified by %. j = j\ + J2 is the total 
electronic angular momentum and J = j + I is the total 
angular momentum of the atom including rotation I. The 
labels m and CI indicate projections along the space-fixed 
Oz axis and intermolecular OZ axis respectively, and w 
represents the gerade (w = 0) and ungerade (w = 1) 
symmetry, respectively, of inversion through the center 
of charge of the molecule. 

The evaluation of the dipole operator must be per- 
formed in the LS basis under the correct symmetry con- 
siderations [|[ , however the inclusion of rji modifies this 
result. In similar fashion to we use the LS basis of 
\jLSCIl,Qs w ), where L = L\ + L2 and S — Si + S2, 
symmetrized according to 



As 



[\rf<yfLSSl L n s )- 



+ (-i) PLS h?7?Lsn L n s }-\ .(i7) 



Note that the coefficients differ from |4( as we include 
spin in the symmetrization. From the above, the matrix 
element of {fji) p d l p in the LS basis is 



{2s2 P \(fji)Pd}\2s2s) 



2V2 



1 



(21) 

where d a t is the atomic dipole moment, w' and Q' L refer 
to the excited state and the condition w — S = even is 
assumed. The presence of (— l) p is a direct result of the 
inclusion of the {fji) p term. 

To convert to the original basis (|16p . the transforma- 
tion (A. 11) of [|| is used. The matrix element of (r)i) p d\ 
in the basis (|16l) reduces to 



2\/2 



1 + (_l) w '+ 1 +i+fl (22) 



where 



F 



Lsn L n s 



^(2S + 1)(2L + l)(2ji + l)(2j 2 + 1) 



x C LSj 

iriLmsmj 




(23) 



{. . .} is a Wigner 9-j coefficient, and dashed quantities 
refer to the 2s2p state and undashed quantities to the 
2s2s state. 
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The matrix element of the rotational part of V™ 1 is 

rot I T~\F* L/,rot\ 



V(2J+1)(2J' + 1) r 



-in 



Using the properties of the rotation matrices, the integral 
can be evaluated [|[ so that 

(vn^AioM,o)io 



2 j/ T j^mjAm^fij/Jff-- l zo 7 



The complete matrix element is therefore 



V mt = - 
eg 



V 

s, W 1 \X r <lpF r ,lpF r< JFJ' n .lF.T 

F/3 



I _|_ (_^iu' + l+J+P 



(26) 



where / is the laser intensity. 

The dipole approximation corresponds to setting kR = 
0, and therefore, since j p {0) = <5 Pl o, to p = in equa- 
tion dm @. The p = term of ' (JUD differs from the 
dipole approximation result by the presence of jo(fci?/2) 
which varies little over the region of the excited bound 
levels but does oscillate and decay significantly outside 
this region. Larger values of p do not contribute strongly 
in the molecular region due to the small magnitude of 
j p (kR/2) for p ^ but are significant in the asymptotic 
region. Lastly, we note that the restriction Aw = ±1 
arising from the dipole approximation is broken, as even 
values of p contribute to couplings to gerade metastablc 
states and odd values of p couple to ungcrade metastablc 
states. 



III. RESULTS 

We first perform the calculations using the same ba- 
sis states for the er~ coupling set that are used in i.e. 
we calculate photoassociation from the gerade metastable 
J = 2 channels to the 0+ , J = 1 adiabatic potential that 
asymptotes to j = 0. Unlike the calculations of Q, the 
asymptotic radiation coupling is zero because the sum- 
mation over p only contributes for < p < 4 and over this 
finite sum the spherical Bcsscl functions j p (kR) — > for 
kR ^> p and kR ^> 1. This means that the dressed state 
formalism or the i?-dependent coupling method need not 
be employed to solve the equations. However, for con- 
sistency between the two calculations, the i?-dependent 
tail-off is introduced exactly as it was used in 0] . 

The photoassociation profiles represent the sponta- 
neous loss from the excited state and are determined 
from the loss of unitarity of the ^-matrix. The photon 



50 



o 40 

i— I 
M— 

o 

ui 30 
'c 

— 20 



u 

01 
l/l 

in 

O 
i— 

U 



10 



0.2 0.4 0.6 0.8 0.2 0.4 

Laser detuning AE (GHz) 



0.6 0.8 



FIG. 1: The photoassociation profiles for intensities of 
64 mW/cm 2 (solid line) and 640 mW/cm 2 (dotted line). The 
left hand plot includes only the p = term of the interac- 
tion coupling and hence only couplings to gerade metastable 
states. The right hand plot includes all terms of the p summa- 
tion and couplings to both gerade and ungerade metastable 
states. 



emission cross section for atoms colliding in the entrance 
channel \a) is 



photon 



A-2 



I'S'a'a I 



(27) 



where a and a' enumerate the open channels and k a 
is the wavenumber of channel a. The photoassociation 
profile is then the average 



photon 



1 



^ ^ ^photon 



(28) 



where n a is the number of open channels. 

Profiles calculated using only the p = contributions 
are shown in the left hand plot of Fig. [TJ The laser in- 
teraction term in this case differs from that in the dipole 
approximation by the presence of jo(kR/2). From the 
plot it is evident that, at large laser intensities, there 
are severe unphysical oscillations in the profiles. Inclu- 
sion of all the terms in the p-summation that contribute 
(p = 0, 2, 4) does not noticeably affect the profiles. 

To mitigate this problem ungerade S = 1 metastable 
states, normally uncoupled in the dipole approximation, 
are introduced as open channels. The p = 1,3 terms of 
the laser coupling contribute to these channels and the 
total loss profiles for this combination of states are pre- 
sented in the right hand plot of Fig. [T] These profiles 
exhibit more expected behavior and possess many simi- 
larities to the original PA profiles found using the dipole 
approximation. 

Closer examination however (see Fig. ^ reveals that 
there is a large difference in behavior between the non- 
dipole and dipole approximation profiles, especially for 
high laser intensities. At low intensities the non-dipolc 
profiles show a reduction in the background loss but the 
resonance peaks have the same shift and width, albeit 
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FIG. 2: The photoassociation profiles for intensities of 
64 mW/cm 2 (solid line), 640 mW/cm 2 (dotted line) and 
3.9 W/cm 2 (dash-dot line). The left hand plot shows profiles 
that result from the use of the dipole approximation whereas 
the right hand plot includes the full summation of the non- 
dipole interaction coupling. 

with an increased resonance strength. At higher laser 
intensities, the resonance parameters still remain identi- 
cal but the resonances in the non-dipole profiles remain 
discernable for larger intensities than those in the dipole 
approximation profiles. In fact, the dipole approxima- 
tion profiles demonstrate a saturation effect as the back- 
ground overwhelms the resonance peaks. Instead, in the 
non-dipole profiles, we observe that the resonance peaks 
decrease in strength such that the overall loss in the 
resonance region is smaller at higher intensities. This 



suggests possible optical suppression of photoassociativc 
processes at high laser intensities. Why this occurs is un- 
known, although it is likely due to destructive interfer- 
ence between the ungerade and gerade metastable chan- 
nels. 

The non-dipole coupling also introduces the possibil- 
ity of more potential photoassociation channels, as the 
restriction | J — 1 1 < J' < \ J + 1 1 is no longer enforced for 
summation terms p > 1. The likelihood that these cou- 
plings induce resonances in the excited state potentials 
is severely suppressed by the spherical Bessel functions 
as j p (kR/2) -C 1 for p ^ and kR -c 1 (i.e. within the 
region of bound levels). This prediction was tested by 
calculating profiles for photoassociation from the J = 2 
metastable levels to the 0+, J = 5 excited levels and, as 
expected, no resonances were visible at any laser inten- 
sity. 

In summary, the PA profiles presented here include 
all possible coupled metastable levels to the 0+, J = 1 
excited state and have been computed without assuming 
the dipole approximation for the laser coupling. The new 
calculations using the non-dipole coupling do not modify 
the resonance parameters that are tabulated in table II of 
Q to the accuracy of the calculations. The observed be- 
havior of the background, and the influence of the dipole 
approximation, are not special to metastable helium but 
are a result of the weakly bound nature of the excited 
levels. It would be of much interest if a similar behav- 
ior were observed in photoassociation of another atomic 
species. 



[1] H. R. Thorsheim, J. Weiner, and P. S. Julienne, Phys. 

Rev. Lett. 58, 2420 (1987) 
[2] J. Weiner, V. S. Bagnato, S. Zilio, and P. S. Julienne, Rev. 

Mod. Phys. 71, 1 (1999) 
[3] D. G. Cocks and I. B. Whittingham, Phys Rev. A 80, 

023417 (2009) 

[4] J. Burke, Ph.D. thesis, University of Colorado, 1999. 



[5] D. M. Brink and G. R. Satchler, Angular Momentum, 2nd 

ed. (Clarendon Press, Oxford, 1968) 
[6] E. E. Nikitin and S. Ya. Umanskii, Theory of Slow Atomic 

Collisions, 1st ed. (Springer, Berlin, 1984) 
[7] A factor of — i is missing from equation (B16) of [3] 



